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We analyze the possibility of probing non-standard neutrino interactions (NSI, for short) through 
the detection of neutrinos produced in a future galactic supernova (SN). We consider the effect of 
NSI on the neutrino propagation through the SN envelope within a three-neutrino framework, paying 
special attention to the inclusion of NSI-induced resonant conversions, which may take place in the 
most deleptonised inner layers. We study the possibility of detecting NSI effects in a Megaton water 
Cherenkov detector, either through modulation effects in the i? e spectrum due to (i) the passage 
of shock waves through the SN envelope, (ii) the time dependence of the electron fraction and (iii) 
the Earth matter effects; or, finally, through the possible detectability of the neutronization v e 
burst. We find that the t> e spectrum can exhibit dramatic features due to the internal NSI-induced 
resonant conversion. This occurs for non-universal NSI strengths of a few %, and for very small 
flavor-changing NSI above a fewxlO . 
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I. INTRODUCTION 

The very first data of the KamLAND collaboration [1| 
have been enough to isolate neutrino oscillations as the 
correct mechanism explaining the solar neutrino prob- 
lem 0, Q] , indicating also that large mixing angle (LMA) 
was the right solution. The 766.3 ton-yr KamLAND data 
sample further strengthens the validity of the LMA os- 
cillation interpretation of the data Q]. 

Current data imply that neutrino have mass. For an 
updated review of the current status of neutrino oscil- 
lations see Q| . Theories of neutrino mass 0, Q typ- 
ically require that neutrinos have non-standard prop- 
erties such as neutrino electromagnetic transition mo- 
ments or non-standard four-Fermi interactions 
(NSI, for short) 0, 0, Q ■ 

The expected magnitude of 
the NSI effects is rather model-dependent. 

Seesaw-type models lead to a non-trivial structure of 
the lepton mixing matrix characterizing the charged and 
neutral current weak interactions The NSI which 
are induced by the charged and neutral current gauge 
interactions may be sizeable 0,0,0,0,0. Alter- 
natively, non-standard neutrino interactions may arise in 
models where neutrinos masses are radiatively "calcula- 
ble" 0, EJ 

Finally, in some supersymmetric unified 
models, the strength of non-standard neutrino interac- 



tions may arise from renormalization and/or threshold 
effects 21 



We stress that non-standard interactions strengths are 
highly model-dependent. In some models NSI strengths 
are too small to be relevant for neutrino propagation, 
because they are either suppressed by some large mass 
scale or restricted by limits on neutrino masses, or both. 
However, this need not be the case, and there are many 
theoretically attractive scenarios where moderately large 
NSI strengths are possible and consistent with the small- 
ness of neutrino masses. In fact one can show that 
NSI may exist even in the limit of massless neutri- 
nos 

0,0,1 16l . Il7l Il8| . Such may also occur in the 



context of fully unified models like SO(10) 22|. 

We argue that, in addition to the precision determi- 
nation of the oscillation parameters, it is necessary to 
test for sub-leading non-oscillation effects that could arise 
from non-standard neutrino interactions. These are nat- 
ural outcome of many neutrino mass models and can be of 
two types: flavor-changing (FC) and non- universal (NU). 
These are constrained by existing experiments (see be- 
low) and, with neutrino experiments now entering a pre- 
cision phase 23j , an improved determination of neutrino 
parameters and their theoretical impact constitute an im- 
portant goal in astroparticle and high energy physics 0| . 

Here we concentrate on the impact of non-standard 
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neutrino interactions on supernova physics. We show 
how complementary information on the NSI parame- 
ters could be inferred from the detection of core-collapse 
supernova neutrinos. The motivation for the study is 
twofold. First, if a future SN event takes place in our 
Galaxy the number of neutrino events expected in the 
current or planned neutrino detectors would be enor- 
mous, O(10 4 — 10 5 ) 24|. Moreover, the extreme con- 
ditions under which neutrinos have to travel since they 
arc created in the SN core, in strongly dclcptonised re- 
gions at nuclear densities, until they reach the Earth, 
lead to strong matter effects. In particular the effect of 
small values of the NSI parameters can be dramatically 
enhanced, possibly leading to observable consequences. 

This paper is planned as follows. In Sec. [IT] we summa- 
rize the current observational bounds on the parameters 
describing the NSI, including previous works on NSI in 
SNe. In Sec. IIIII we describe the neutrino propagation 
formalism as well as the SN profiles which will be used. 
In Sec. lIVI we analyze the effect of NSI on the v propaga- 
tion in the inner regions near the neutrinosphere and in 
the outer regions of the SN envelope. In Scc.|V]we discuss 
the possibility of using various observables to probe the 
presence of NSI in the neutrino signal of a future galactic 
SN. Finally in Sec. I VII we present our conclusions. 
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32j. On the other hand the analysis 



-> vv^ cross section measured at LEP II 
leads to a bound on < 0.5 [33j. Future prospects 

to improve the current limits imply the measurement of 
sin 2 dw leptonically in the scattering off electrons in the 
target, as well as in neutrino deep inelastic scattering in 
a future neutrino factory. The main improvement would 
be in the case of |e£f| and |e£T|, where values as small 
as 10~ 3 and 0.02, respectively, could be reached [31 1. 



The search for flavor violating processes involving 
charged leptons is expected to restrict corresponding neu- 
trino interactions, to the extent that the SU(2) gauge 
symmetry is assumed. However, this can at most give 
indicative order-of-magnitude restrictions, since we know 
SU(2) is not a good symmetry of nature. Using radiative 
corrections it has been argued that, for example, /i — e 
conversion on nuclei like in the case of \i~Ti also con- 



31] 



strains \ef e \ < 7.7 x 10~ 4 

Non-standard interactions can also affect neutrino 
propagation through matter, probed in current neutrino 
oscillation experiments. The bounds so obtained apply to 
the vector coupling constant of the NSI, s^X = ^Jn+^Jg, 
since only this appears in neutrino propagation in mat- 
ter [91L 



II. PRELIMINARIES 



B. Solar and reactor 



A large class of non-standard interactions may be 
parametrized with the effective low-energy four-fermion 
operator: 



C 



NSI 



%2y/2G F (v a <y li Lu l3 )(J>fPf) , (1) 



where P = L, R and / is a first generation fermion: 
e, u, d. The coefficients denote the strength of the 
NSI between the neutrinos of flavors a and f3 and the 
P— handed component of the fermion /. 

Current constraints on e a o come from a variety of dif- 
ferent sources, which we now briefly list. 

A. Laboratory 



Neutrino scattering experiments 
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29[ provide the following bounds, \slE\ < 10 



The role of neutrino NSIs as sublcading effects on the 
solar neutrino oscillations and KamLAND has been re- 
cently considered in Rcf. 
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with the following 
bounds at 90 % CL for e = — sin-i^^er with the al- 
lowed range —0.93 < e < 0.30, while for the diagonal 
term e' = sin 2 $ 2 3£tt — £ fe \ the only forbidden region is 



[0.20,0.78] [36|. Only in the ideal case of infinitely pre- 
cise solar neutrino oscillation parameters determination, 
the allowed range would "close from the left" for negative 
NSI parameter values, at —0.6 for e and —0.7 for e'. 

C. Atmospheric and accelerator neutrinos 

Non-standard interactions involving muon neutrinos 
can be constrained by atmospheric neutrino experiments 
as well as accelerator neutrino oscillation searches at 
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K2K and MINOS. In Rcf. [37| Supcr-Kamiokande and 
MACRO observations of atmospheric neutrinos were con- 
sidered in the framework of two neutrinos. The limits ob- 
tained were -0.05 < ef T < 0.04 and - e*%\ < 0.17 
at 99 % CL. The same data set together with K2K were 
recently considered in Refs. 



MX 



to study the nonstan- 
dard neutrino interactions in a three generation scheme 
under the assumption e efl = £ w = e^ T = 0. The al- 
lowed region of e TT obtained for values of e er smaller 



than ©(lO- 1 ) becomes Y, f = uAe e^N f /N e < 0.2 [39| , 
where Nf stands for the fermion number density. 

D. Cosmology 

If non-standard interactions with electrons were large 
they might also lead to important cosmological and as- 
trophysical implications. For instance, neutrinos could 
be kept in thermal contact with electrons and positrons 
longer than in the standard case, hence they would share 
a larger fraction of the entropy release from e annihi- 
lations. This would affect the predicted features of the 
cosmic background of neutrinos. As recently pointed out 
in Ref 40( required couplings are, though, larger than 



the current laboratory bounds. 



E. NSI in Supernovae 

According to the currently accepted supernova (SN) 
paradigm, neutrinos are expected to play a crucial role 
in SN dynamics. As a result, SN physics provides 
a laboratory to probe neutrino properties. Moreover, 
many future large neutrino detectors are currently be- 
ing discussed 4l|. The enormous number of events, 
O(10 4 — 10 5 ) that would be "seen" in these detectors in- 
dicates that a future SN in our Galaxy would provide a 
very sensitive probe of non-standard neutrino interaction 
effects. 

The presence of NSI can lead to important conse- 
quences for the SN neutrino physics both in the highly 
dense core as well as in the envelope where neutrinos 
basically freely stream. 

The role of non-forward neutrino scattering processes 
on heavy nuclei and free nucleons giving rise to flavor 
change within the SN core has been recently analyzed in 



Ref. 42j, |43j . The main effect found was a reduction in 
the core electron fraction Y e during core collapse. A lower 
Y e would lead to a lower homologous core mass, a lower 
shock energy, and a greater nuclear photon-disintegration 
burden for the shock wave. By allowing a maximum 
AY e = —0.02 it has been claimed that e ea < 10 -3 , where 
a = n, t 0' 

On the other hand it has been noted since long ago 
that the existence of NSI plays an important role in the 
propagation of SN neutrinos through the envelope lead- 
ing to the possibility of a new resonant conversion. In 



44, |45j it would 
131. Two basic 



contrast to the well known MSW effect 
take place even for massless neutrinos 
ingredients are necessary: universal and flavor changing 
NSI. In the original scheme neutrinos were mixed in the 
leptonic charged current and universality was violated 
thanks to the effect of mixing with heavy gauge singlet 
leptons 0, [ljj . Such resonance would induce strong neu- 
trino flavor conversion both for neutrinos and antineutri- 
nos simultaneously, possibly affecting the neutrino sig- 
nal of the SN1987A as well as the possibility of having 
r— process nucleosynthesis. This was first quantitatively 
considered within a two- flavor v e — v r scheme, and bounds 
on the relevant NSI parameters were obtained using both 
arguments 
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One of the main features of the such "internal" or 
"massless" resonant conversion mechanism is that it re- 
quires the violation of universality, its position being 
determined only by the matter chemical composition, 
namely the value of the electron fraction Y ei and not by 
the density. In view of the experimental upper bounds 
on the NSI parameters such new resonance can only take 
place in the inner layers of the supernova, near the neu- 
trinosphere, where Y e takes its minimum values. In this 
region the values of Y e are small enough to allow for 
resonance conversions to take place in agreement with 
existing bounds on the strengths of non-universal NSI 
parameters. 

The SN physics implications of another type of NSI 
present in supcrsymmctric R- par ity violating models 
have also been studied in Ref. 47J, again for a system 



of two neutrinos. For definitcness NSI on d— quarks were 
considered, in two cases: (i) massless neutrinos without 
mixing in the presence of flavor-changing (FC) and non- 
universal (NU) NSIs, and (ii) neutrinos with eV masses 
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and FC NSI. Different arguments have been used in 
order to constrain the parameters describing the NSI, 
namely, the SN1987A signal, the possibility to get suc- 
cessful r— process nucleosynthesis, and the possible en- 
hancement of the energy deposition behind the shock 
wave to reactivate it. 



On the other hand several subsequent articles 48|, |49|, 
[5q| | considered the effects of NSI on the neutrino propa- 
gation in a three-neutrino mixing scenario for the case 
Y e > 0.4, typical for the outer SN envelope. Together 
with the assumption that e^i < 10~ 2 this prevents the 
appearance of internal resonances in contrast to previous 
references. 

Motivated by supersymmetric theories without R par- 



ity, in Ref. 



the authors considered the effects of 



small-strength NSI with d— quarks. Following the for 



malism developed in Rcfs. 



51 



they studied the cor- 



rections that such NSI would have on the expressions 
for the survival probabilities in the standard resonances 
MSW-H and MSW-L. A similar analysis was performed 
in Rcf. [49! assuming Z-induced NSI interactions orig- 
inated by additional heavy neutrinos. A phenomcno- 
logic al g eneralization of these results was carried out in 
Ref. [501 ] . The authors found an analytical compact ex- 
pression for the survival probabilities in which the main 
effects of the NSI can be embedded through shifts of the 
mixing angles #x2 an d $13. In contrast to similar expres- 
sions found previously these directly apply to all mixing 
angles, and in the case with Earth matter effects. The 
main phenomenological consequence was the identifica- 
tion of a degeneracy between $13 and e ea , similar to the 
analogous "confusion" between $13 and the correspond- 
ing NSI parameter noted to exist in the context of long- 
baseline neutrino oscillations 



r,:-i 



54]. 



We have now re-considered the general three-neutrino 
mixing scenario with NSI. In contrast to previous 



work 



48 



49 



0, 



we have not restricted ourselves to 



large values of Y e , discussing also small values present 
in the inner layers. This way our generalized descrip- 
tion includes both the possibility of neutrinos having the 
"massless" NSI-induced resonant conversions in the in- 
ner layers of the SN envelope 
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4G 



as well as the 



"outer" oscillation-induced conversions [48 



49 



50| [92|. 



III. NEUTRINO EVOLUTION 

In this section we describe the main ingredients of our 
analysis. Our emphasis will be on the use of astrophys- 
ically realistic SN matter and Y e profiles, characterizing 
its density and the matter composition. Their details, 
in particular their time dependence, are crucial in deter- 
mining the way the non-standard neutrino interactions 
affect the propagation of neutrinos in the SN medium. 

A. Evolution Equation 

As discussed in Sec. [II] in an unpolarizcd medium the 
neutrino propagation in matter will be affected by the 
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vector coupling constant of the NSI, e~X 
The way the neutral current NSI modifies the neu- 
trino evolution will be parametrized phenomenologically 
through the effective low-energy four-fermion operator 
described in Eq. ([T]). We also assume e^g S 3?, neglect- 
ing possible CP violation in the new interactions. 

Under these assumptions the Hamiltonian describing 
the SN neutrino evolution in the presence of NSI can be 
cast in the following form 94 1 



l-T-Va = (H kin 
or 



where -f/kin stands for the kinetic term 



2E 



(2) 



(3) 



with M 2 = diag(mf, rrio, TO3), and U the three-neutrino 
lepton mixing matrix [6( in the PDG convention 55j and 
with no CP phases. 

The second term of the Hamiltonian accounts for the 
interaction of neutrinos with matter and can be split into 
two pieces, 



H 



int 



std 



Hi 



(4) 



The first term, H^f describes the standard interaction 
with matter and can be written as Hf*£ = diag (Vcc, 0, 0) 
up to one loop corrections due to different masses of the 
muon and tau leptons 56J. The standard matter poten- 



tial for neutrinos is given by 

Vcc = V2G F N e 



V oP Y e 



(5) 
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where Vq w 7.6 x 10~ 14 cV, the density is given in g/cm 3 , 
and Y e stands for the relative number of electrons with 
respect to baryons. For antineutrinos the potential is 
identical but with the sign changed. 

The term in the Hamiltonian describing the non- 
standard neutrino interactions with a fermion / can be 
expressed as, 

(h™u = £ <yL)<#, (6) 

f— e,u,d 

with (V^ si ) a p = VZGpNfe-^p. For definiteness and mo- 
tivated by actual models, for example, those with broken 
R parity supersymmetry we take for / the down-type 
quark. However, an analogous treatment would apply 
to the case of NSI on up-type quarks, the existence of 
NSI with electrons brings no drastic qualitative differ- 
ences with respect to the pure oscillation case (see be- 
low). Therefore the NSI potential can be expressed as 
follows, 



(V n d si U = e d aP V p(2 - Y e ) 



(7) 



From now on we will not explicitely write the superindcx 
d. In order to further simplify the problem we will rede- 



fine the diagonal NSI parameters so that 



0. as one 



can easily see that subtracting a matrix proportional to 
the identity leaves the physics involved in the neutrino 
oscillation unaffected. 

B. Supernova matter profiles 

Neutrino propagation depends on the supernova mat- 
ter and chemical profile through the effective potential. 
This profile exhibits an important time dependence dur- 
ing the explosion. Fig.[T]shows the density p(t, r) and the 
electron fraction Y e (t, r) profiles for the SN progenitor as 
well as at different times post-bounce. 

Progenitor density profiles can be roughly 
parametrized by a power-law function 



(8) 



where po ~ 10 4 g/cm 3 , Rq ~ 10 9 cm, and n ~ 3. The 
electron fraction profile varies depending on the matter 
composition of the different layers. For instance, typical 
values of Y e between 0.42 and 0.45 in the inner regions 



are found in stellar evolution simulations 57 | . In the in- 
termediate regions, where the MSW H and L-rcsonanccs 
take place Y e sa 0.5. This value can further increase in 
the most outer layers of the SN envelope due to the pres- 
ence of hydrogen. 

After the SN core bounce the matter profile is affected 
in several ways. First note that a front shock wave starts 
to propagate outwards and eventually ejects the SN enve- 
lope. The evolution of the shock wave will strongly mod- 
ify the density profile and therefore the neutrino propa- 
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gation 58|, |59| . Following Ref. |60j we shall assume that 
the structure of the shock wave is more complicated and 
an additional "reverse wave" appears due to the collision 
of the neutrino-driven wind and the slowly moving mate- 
rial behind the forward shock, as seen in the upper panel 
of Fig. [rj 0- 

On the other hand, the electron fraction is also affected 
by the time evolution as the SN explosion proceeds. Once 
the collapse starts the core density grows so that the neu- 
trinos become eventually effectively trapped within the 
so called "neutrinosphere" . At this point the trapped 
electron fraction has decreased until values of the order of 
0.33 [6lj] . When the inner core reaches the nuclear density 
it can not contract any further and bounces. As a con- 
sequence a shock wave forms in the inner core and starts 
propagating outwards. When the newly formed super- 
nova shock reaches densities low enough for the initially 
trapped neutrinos to begin streaming faster than the 
shock propagates 62|, a breakout pulse of v e is launched. 
In the shock-heated matter, which is still rich of elec- 
trons and completely disintegrated into free neutrons and 
protons, a large number of v e are rapidly produced by 
electron captures on protons. They follow the shock on 
its way out until they are released in a very luminous 
flash, the breakout burst, at about the moment when 
the shock penetrates the neutrinosphere and the neutri- 
nos can escape essentially unhindered. As a consequence, 
the lepton number in the layer around the neutrinosphere 
decreases strongly and the matter neutronizes Q. The 
value of Y e steadily decreases in these layers until val- 
ues of the order of 0(1O -2 ). Outside the neutrinosphere 
there is a steep rise until Y e m 0.5. This is a robust 
feature of the neutrino-driven baryonic wind. Neutrino 
heating drives the wind mass loss and causes Y e to rise 
within a few 10 km from low to high values, between 0.45 
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and 0.55 64(, see bottom panel of Fig. [TJ Inspired in the 



numerical results of Ref. 60( we have parametrized the 
behavior of the electron fraction near the ncutrinosphcrc 
phcnomcnologically as, 



Y e = a + 6arctan[(r — ro)/r s ] , 



(9) 



where a « 0.23 - 0.26 and b w 0.16 - 0.20. The param- 
eters ro and r s describe where the rise takes place and 
how steep it is, respectively. As can be seen in Fig. [1] 
both decrease with time. 




FIG. 1: Density (upper panel) and electron fraction (bottom 
panel) profiles for the SN progenitor and at different instants 
after the core bounce, from Ref. The regions where the 

H (yellow) and the L (cyan) resonance take place are also 
indicated, as well as the NSI-induced / (gray) resonance for 
the parameters e ee = 0, e T t < 0.07 and |e MT | < 0.05 



IV. THE TWO REGIMES 

In order to study the neutrino propagation through 
the SN envelope we will split the problem into two differ- 
ent regions: the inner envelope, defined by the condition 



2 

atm 



777-2, an d the 



Vcc > Aml tm /(2E) with Am 
outer one, where Amj tm /(2£) > Vcc- From the upper 
panel of Fig. [1] one can see how the boundary roughly 



varies between r w 10 s cm and 10 9 cm, depending on the 
time considered. This way one can fully characterize all 
resonances that can take place in the propagation of su- 
pernova neutrinos, both the outer resonant conversions 
related to neutrino masses and indicated as the upper 
bands in Fig. [TJ and the inner resonances that follow 
from the presence of non-standard neutrino interactions, 
indicated by the band at the bottom of the same figure. 
Here we pay special attention to the use of realistic mat- 
ter and chemical supernova profiles and three-neutrino 
flavors thus generalising previous studies. 

A. Neutrino Evolution in the Inner Regions 

Let us first write the Hamiltonian in the inner layers, 
where Hi nt 3> i?km- In this case the Hamiltonian can be 
written as 



V p(2 - Y e ) 



I 



\ 



2-Y e ^ tee 













£ er 


£{1T 





(10) 

When the value of the e a p is of the same order as the 
electron fraction Y e internal resonances can arise 
Taking into account the current constraints on the e's 
discussed Jn Sec. |H] one sees that small values of Y e are 
required 
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471 ]. As a result, these can only take place 
in the most deleptonised inner layers, close to the neu- 
trinosphere, where the kinetic terms of the Hamiltonian 
are negligible. 

Given the large number of free parameters e a p in- 
volved we consider one particular case where |e eAl | and 
\e eT \ are small enough to neglect a possible initial mixing 
between v e and or v T . Barring fine tuning, this basi- 
cally amounts to |e e ^|, |£er| ^ 10 -2 . According to the 
discussion of Sec.|n]ee/i automatically satisfies the condi- 
tion, whereas one expects that the window |e eT | > 10~ 2 
will eventually be probed in future experiments. 

Since the initial fluxes of and v T arc expected to be 
basically identical, it is convenient to redefine the weak 
basis by performing a rotation in the [i — t sector: 

/ v e 
v' 




S23' C 2 3' 



V 



111) 
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where C23' and S23' stand for cos(t9 23 ) and sin(i9 23 ), re- 
spectively. The angle z? 23 can be written as 



tan(2tf 2 



2H- 



2.1 



23 / 



133 £tt 

The Hamiltonian becomes in the new basis 



H 'af3 



V p(2 - Y t 



1 2-Y- 



(12) 

(13) 
,(14) 



V 



where 



- e// 



£e/xC23' 



-S23' 



^ en S23> + £erC23' 



{e TT -Je% T + Ael T )/2 



= ( £tt + ,/4 t + 44 t )/2 



(15) 
(16) 
(17) 

(18) 



With our initial assumptions on e ea one notices that 
the new basis v' a basically diagonalizes the Hamiltonian, 
and therefore coincides roughly with the matter eigen- 
state basis. A novel resonance can arise if the condition 
H' ee = H' TT is satisfied, we call this /-resonance, / stand- 
ing for "internal" [3] . The corresponding resonance con- 
dition can be written as 



Y 1 



2s 1 



1 



(19) 



where e 1 is defined as e' TT — e ee . In Fig. [5] we represent 
the range of e ee and e' TT leading to the /-resonance for 
an electron fraction profile between different Y™ ln, s and 
y-max _ g g j^. j g i m p 0r tant to notice that the value 
of y e mm depends on time. Right before the collapse the 
minimum value of the electron fraction is around 0.4. 
Hence the window of NSI parameters that would lead 
to a resonance would be relatively narrow, as indicated 
by the shaded (yellow) band in Fig. O As time goes on 
y e min decreases to values of the order of a few %, and 
as a result the region of parameters giving rise to the /- 
resonance significantly widens. For example, in the range 
\e ee \ < 10 -3 possibly accessible to future experiments one 
sees that the /-resonance can take place for values of e' TT 
of the order of O(10 -2 ). This indicates that the potential 
sensitivity on NSI parameters that can be achieved in su- 
pernova studies is better than that of the current limits. 
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FIG. 2: Contours of Y e as function of e ee and e' TT accord- 
ing to Eq. (|19|l for different values of Y e . The region in yel- 
low represents the region of parameters that gives rise to I- 
resonance before the collapse. The arrows indicate how this 
region widens with time. 



As seen in Fig. Q] in order to fulfill the /-resonance con- 
dition for such small values of the NSI parameters the 
values of Y e must indeed lie, as already stated, in the 
inner layers. 

Several comments are in order: First, in contrast to 
the standard H and /-resonances, related to the kinetic 
term, the density itself does not explicitly enter into the 
resonance condition, provided that the density is high 
enough to neglect the kinetic terms. Analogously the en- 
ergy plays no role in the resonance condition, which is 
determined only by the electron fraction Y e . Moreover, 
in contrast to the standard resonances, the /-resonance 
occurs for both neutrinos and antineutrinos simultane- 
ously (l3j [. Finally, as indicated in Fig. [3]the v^s (D e ) are 
not created as the heaviest (lightest) state but as the in- 
termediate state, therefore the flavor composition of the 
neutrinos arriving at the //-resonance is exactly the op- 
posite of the case without NSI. As we show in Scc.|Vl this 
fact can lead to important observational consequences. 

In order to calculate the hopping probability between 
matter eigenstates at the /-resonance we use the Landau- 
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Zencr approximation for two flavors 



FIG. 3: Level-crossing schemes, first panel is for the case of 
normal hierarchy (oscillations only), the second includes the 
NSI effect. The two lower panels correspond to the inverse hi- 
erarchy, oscillations only and oscillations + NSI, respectively. 



P 



LZ 



(20) 



where jj stands for the adiabaticity parameter, which 
can be generally written as 



Ef - Ef 



2tf n 



(21) 



where i? m = d^ m /dr. If one applies this for- 
mula to the e — t' box of Eq. (fT4|) assuming that 



tan2i?7 = 2H' eT /{H' TT - H ee ) 

1 1/2 , 

one gets 



and E r 2 n - El 



[{H' TT -H ee f 



AH' 2 









ri 


(1 + e'YYe 


ri 



(H' TT - H e . 
4 x 10 9 r- B pne* /(e 1 ) . 



(22) 



where the parametrization of the Y e profile has been de- 
fined as in Eq. ((9]) with b = 0.16. The density pu rep- 
resents the density in units of 10 11 g/cm 3 , r Sj 5 stands 
for r s in units of 10 5 cm, and fie 1 ) is a function whose 
value is of the order 0(1) in the range of parameters we 
are interested in. Taking all these factors into account 
it follows that the internal resonance will be adiabatic 
provided that e' eT > 10~ 5 , well below the current limits, 
in full numerical agreement with, e. g., Ref. 47 1. 

In Fig. [4] we show the resonance condition as well as 
the adiabaticity in terms of e TT and e eT assuming the 
other e a fj = 0. In order to illustrate the dependence on 
time we consider profiles inspired in the numerical profiles 
of Fig. Q] at t = 2 s (upper panel) and 15.7 s (bottom 
panel). For definiteness we take y^ min as the electron 
fraction at which the density has value of 5 x 10 11 g/cm 3 . 
For comparison with Fig. [2] we have assumed Y e mm = 
10~ 2 in the case of 15.7 s. We observe how the border 
of adiabaticity depends on e TT through the value of the 
density at rj which in turn depends on time. 

Before moving to the discussion of the outer resonances 
a comment is in order, namely, how does the formalism 
change for other non-standard interaction models. First 
note that the whole treatment presented above also ap- 
plies to the case of NSI on up- type quarks, except that 
the position of the internal resonance shifts with respect 
to the down-quark case. Indeed, in this case the NSI 
potential 



(Ci) 



a/3 



-a/3 



V p{l + Y e ) 



(23) 
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FIG. 4: Contours of constant jump probability at the I- 
resonance in terms of e TT and e eT for two profiles correspond- 
ing to Fig.[T]at 2 s with a — 0.235 and b — 0.175 (upper panel) 
and 15.7 s with a = 0.26 and b = 0.195 (bottom panel). For 
simplicity the other e's have been set to zero. 

would induce a similar internal resonance for the condi- 
tion Y e = e' /{l-e 1 ). 

In contrast, for the case of NSI with electrons, the 
NSI potential is proportional to the electron fraction, and 
therefore no internal resonance would appear. 



B. Neutrino Evolution in the Outer Regions 

In the outer layers of the SN envelope neutrinos can un- 
dergo important flavor transitions at those points where 
the matter induced potential equals the kinetic terms. 
In absence of NSI this condition can be expressed as 
Vcc ~ Am 2 /(2£). Neutrino oscillation experiments in- 
dicate two mass scales, ATO 2 tm and AttIq = m 2 — m\ 
hence two different resonance layers arise, the so-called 
//-resonance and the L-resonance, respectively. 

The presence of NSI with values of \e a i3\ < 10 -2 modi- 



fies the properties of the H and L transitions 48|, [49J, [5C | . 
In particular one finds that the effects of the NSI can be 
described as in the standard case by embedding the e's 
into effective mixing angles [5fJ. An analogous "confu- 
sion" between sin §i 3 and the corresponding NSI param- 
eter e eT has been pointed out in the context of long- 
baseline neutrino oscillations in Rcfs. 53 



In this section we perform a more general and com- 
plementary study for slightly higher values of the NSI 
parameters: \e a p\ > few 1CT 2 , still allowed by current 
limits, and for which the /-resonance could occur. 

The phenomenological assumption of hierarchical 
squared mass differences, |Am 2 tm | ^> Am|,, allows, for 
not too large e's, a factorization of the dynamics 
into two 2v subsystems roughly decoupled for the H 
and L transitions 65[. To isolate the dynamics of the 



H transition, one usually rotates the neutrino flavor ba- 
sis by UH'&2?,\ and extracts the submatrix with indices 



5fJ. Whereas this method works perfectly for 
small values of e a p it can be dangerous for values above 
1CP 2 . In order to analyze how much our case deviates 
from the simplest approximation we have performed a 
rotation with the angle # 2 3 = ^23 — ol instead of just 
$23- By requiring that the new rotation diagonalizes the 
submatrix (2,3) at the //-resonance layer one obtains the 
following expression for the correction angle a 



tan(2a) = [A Q s2 12 s 13 + V T N T SI s2 



2:5 



(A atm + iA )c 2 3 + iA Q c2 12 (-3 + c2 13 ) 



where A a 



+V T N T SI c2 23 + 2V» SI S 2 23 ] 
= Am 2 tm /(2£) and A e 



(24) 



Am 2 Q /(2E). In 
our notation Sij and s2^ represent sin-i^- and sin(2$jj), 
respectively. The parameters Cy and c2y are analogously 
defined. In the absence of NSI a is just a small correction 
to t?23 



97] 



tan(2a) » A s2 12 s 13 /A atm c 2 3 < O(10" 3 ) 



(25) 



In order to calculate a we need to know the H- 
resonance point. To calculate it one can proceed as in 
the case without NSI, namely, make the $ 23 rotation and 
analyze the submatrix (1,3). The new Hamiltonian H'^o 
has now the form 

Ke = ^0PK+£ee(2-i;)]+AatmS 2 3 
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+ A 0( C ?3 S 12 + s 13) > 

h'; t = K, P (2-y e )^ T + A atm 44 

+ Aq [cl 3 C 2 a + (s Q Ci2 + C a Si 2 Si 3 ) 2 ] , 

H'J T = V p(2-Y e )e'l r + ^A atm s2 13 c a 

+ ^A Q (-c 13 s a s2 12 + c( 2 c a s2 13 ) . (26) 
We have defined e" T = e TT c^ 3 _ a + e^ T s2 23 ^ a , and e" r = 

£eTC 2 3-a+£e^S 23 ^ a , where S 23 _ a = Smfe-a), C23- a = 

cos(# 23 - a), and s2 23 _ Q = sin(2i9 23 - 2a), c2 23 _ Q = 
cos(2^ 2 3 — 2a). The resonance condition for the H tran- 
sition, H " e = iJ" T can be then written as 

V 0P H [Y e H + (e ee - 4' T )(2 - n ff )] = Aat m (c? 3 4 " 40 

+ A [ C 2 2 (c? 3 - C 2 S 2 3 ) _ S 2 S 2 2 + 1 s2qS 2 12Si3 ](27) 

It can be easily checked how in the limit of e a p — > one 
recovers the standard resonance condition, 



V p"Y e 



A a t m c2 



13 



(28) 



In the region where the iJ-resonance occurs Yj 1 w 0.5. 

Taking into account Eqs. (|24|) and (|27|) one can already 
estimate how the value of a changes with the NSI param- 
eters. In Fig. Owe show the dependence of a on the e TT 
after fixing the value of the other NSI parameters. One 
can see how for e TT > 10~ 2 the approximation of neglect- 
ing a significantly worsens. Assuming fl 23 = 7r/4 and a 
fixed value of £ MT one can easily see that e TT basically 
affects the numerator in Eq. ([24]) . Therefore one expects 
a rise of a as the value of e TT increases, as seen in Fig. [5] 
The dependence of a on e^ T is correlated to the rela- 
tive sign of the mass hierarchy and e Mr . For instance, 
for normal mass hierarchy and positive values of e^ T the 
dependence is inverse, namely, higher values of e^ T lead 
to a suppression of a. Apart from this general behav- 
ior, a also depends on the diagonal term e ee as seen in 
Fig. [5] This effect occurs by shifting the resonance point 
through the resonance condition in Eq. (|27[). 

One can now calculate the jump probability be- 
tween matter eigenstates in analogy to the /-resonance 
by means of the Landau-Zencr approximation, sec 
Eqs. 0, (EH), and [H 



H 

LZ 



(29) 
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FIG. 5: Angle a as function of e TT for different values of e ee 
and e MT , in the case of neutrinos of energy 10 MeV, with 
normal mass hierarchy, and sf 3 = 10 -5 . The other NSI pa- 
rameters take the following values: e eM = and e eT = 10~ 3 . 



ii-resonancc, which can be written as 



1h 



4if' 



(30) 



<H" - H" ) 

where the expressions for H'^ are given in Eqs (f2"6"|) . 

Let us first consider the case \e a p\ < 10~ 2 . In this case 
a ?s and one can rewrite the adiabaticity parameter as 

A atm sin 2 (2^^ / ) 



7ff 



cos(2i?5{ / )|dlnV/dr| r 



where 



"13 



d lz + e'i T {2-Y e )/Y e 



(31) 



(32) 



where represents the adiabaticity parameter at the 



in agreement with Ref. 50(. For slightly larger e's there 
can be significant differences. In Fig. [5] we show P^ z in 
the e eT -e TT plane for antineutrinos with energy 10 MeV 
in the case of inverse mass hierarchy, using Eq. (|2"9"|) with 
(upper panel) and without (bottom panel) the a cor- 
rection. The values of z?i 3 and e eT have been chosen so 
that the jump probability lies in the transition regime be- 
tween adiabatic and strongly non adiabatic. In the limit 
of small £ TT , a becomes negligible and therefore both re- 
sults coincide. From Eq. (|3"Tj) one sees how as the value of 
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e eT increases ju gets larger and therefore the transition 
becomes more and more adiabatic. For negative values 
of e eT there can be a cancellation between e eT and z? 



13: 



and as a result the transition becomes non-adiabatic. 



An additional consequence of Eq. is that a degen- 
eracy between e eT and $13 arises. This is seen in Fig. [7j 
which gives the contours of P^z m terms of e eT and $13 
for e TT = 1CP 4 . One sees clearly that the same Landau- 
Zener hopping probability is obtained for different com- 
binations of e eT and $13. This leads to an intrinsic "con- 
fusion" between the mixing angle and the corresponding 
NSI parameter, which can not be disentangled only in 
the context of SN neutrinos, as noted in Ref. 5C|. 



We now turn to the case of \s TT \ > 10~ 2 . As |e TT | 
increases the role of a becomes relevant. Whereas in the 
bottom panel Pjjz remains basically independent of e TT , 
one can see how in the upper panel P^ z becomes strongly 
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One sees that for positive values of e TT it tends to adi- 
abaticity whereas for negative values to non-adiabaticity. 
This follows from the dependence of H'J r on a, essen- 
tially through the term — A q ci3Sq,s2i2, see Eq. (|26|). For 
\e TT \ > 10 -2 one sees that sin a starts being important, 
and as a result this term eventually becomes of the same 
order as the others in H" T . At this point the sign of 
e TT , and so the sign of sin a, is crucial since it may con- 
tribute to the enhancement or reduction of H'J T . This 
directly translates into a trend towards adiabaticity or 
non-adiabaticity, seen in Fig. [6l Thus, for the range of 
e TT relevant for the NSI-induced internal resonance the 
adiabaticity of the outer H resonance can be affected in 
a non-trivial way. 

Turning to the case of the L transition a similar expres- 
sion can be obtained by rotating the original Hamiltonian 



by u(-& 13 yu(-& 23 ) 



48 



3. 



However, in contrast to the 



case of the P-resonance, where the mixing angle $13 is 
still unknown, in the case of the L transition the angle 
$12 has been shown by solar and reactor neutrino exper- 
iments to be large As a result, for the mass scale A Q 
this transition will always be adiabatic irrespective of the 
values of £ a| g, and will affect only neutrinos. 



FIG. 6: Landau-Zener jump probability isocontours at the H- 
resonance in terms of e e r and e TT for 10 MeV antineutrinos 
in the case of inverted mass hierarchy. Upper panel: a given 
by Eq. (|24[) . Bottom panel: a set to zero. The remaining 
parameters take the following 
10 -3 , See = £eu = 0. See text. 



V. OBSERVABLES AND SENSITIVITY 

As mentioned in the introduction one of the major mo- 
tivations to study NSI using the neutrinos emitted in a 
SN is the enhancement of the NSI effects on the neutrino 
propagation through the SN envelope due to the specific 
extreme matter conditions that characterize it. In this 
section we analyze how these effects translate into ob- 
servable effects in the case of a future galactic SN. 

Schematically, the neutrino emission by a SN can be di- 
vided into four stages: Infall phase, neutronization burst, 
accretion, and Kelvin-Helmholtz cooling phase. During 
the infall phase and neutronization burst only v^s are 
emitted, while the bulk of neutrino emission is released 
in all flavors in the last two phases. Whereas the neutrino 
emission characteristics of the two initial stages are basi- 
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FIG. 7: Landau-Zener jump probability isocontours at the 
//-resonance in terms of e eT and $13 for e TT = 1CP 4 . An- 
tineutrinos with energy 10 MeV and inverted mass hierarchy 
has been assumed. 



cally independent of the features of the progenitor, such 
as the core mass or equation of state (EoS), the details 
of the neutrino spectra and luminosity during the ac- 
cretion and cooling phases may significantly change for 
different progenitor models. As a result, a straightfor- 
ward extraction of oscillation parameters from the bulk 
of the SN neutrino signal seems hopeless. Only features 
in the detected neutrino spectra which are independent 
of unknown SN parameters should be used in such an 
analysis 



The question then arises as to how can one obtain in- 
formation about the NSI parameters. Taking into ac- 
count that the main effect of NSI is to generate new in- 
ternal neutrino flavor transitions, one possibility is to in- 
voke theoretical arguments that involve different aspects 
of the SN internal dynamics. 

In Ref. 47j it was argued that such an internal flavor 



conversion during the first second after the core bounce 
might play a positive role in the so-called SN shock re- 
hcating problem. It is observed in numerical simula- 



G7 



G8 



G9 



70j that as the shock wave propa- 



tions 

gates it loses energy until it gets stalled at a few hun- 
dred km. It is currently believed that after neutrinos 
escape the SN core they can to some extent deposit en- 
ergy right behind and help the shock wave continue out- 



wards. On the other hand it is also believed that due to 
the composition in matter of the protoncutronstar (PNS) 
the mean energies of the different neutrino spectra obey 
(E Ve ) < (Ep e ) < (E v v T )- This means that a reso- 
nant conversion between v e {^e) and z^, T (£^, r ) between 
the neutrinosphere and the position of the stalled shock 
wave would make the v e (v e ) spectra harder, and there- 
fore the energy deposition would be larger, giving rise to 
a shock wave regeneration effect. 

Another argument used in the literature was the pos- 
sibility that the r— process nucleosynthesis, responsible 
for synthesizing about half of the heavy elements with 
mass number A > 70 in nature, could occur in the region 
above the neutrinosphere in SNe [ll[ 0| ■ A necessary 
condition is Y e < 0.5 in the nucleosynthesis region. The 
value of the electron fraction depends on the neutrino 
absorption rates, which arc determined in turn by the 
v e {pe) luminosities and energy distribution. These can 
be altered by flavor conversion in the inner layers due to 
the presence of NSI. Therefore by requiring the electron 
fraction be below 0.5 one can get information about the 
values of the NSI parameters. 

While it is commonly accepted that neutrinos will play 
a crucial role in both the shock wave re-heating as well 
as the r— process nucleosynthesis, there are still other as- 
trophysical factors that can affect both. While the issue 
remains under debate we prefer to stick to arguments 
directly related to physical observables in a large water 
Cherenkov detector. There are several possibilities. 



(A) the modulations in the v e spectra due to the 



sag e of shock waves through the supernova 
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(B) the modulation in the v e spectra due to the time 
dependence of the electron fraction, induced by the 
/-resonance 

(C) the modulations in the v e spectra due to the Earth 
matter 



(D) detectability of the neutronization v e burst [77I. [t8| 

Three of these observables, 1, 3 and 4 have already been 
considered in the literature in the context of neutrino 
oscillations. Here we discuss the potential of the above 
promising observables in providing information about the 
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TABLE I: Definition of the neutrino schemes considered in 
terms of the hierarchy, the value of $13, and the presence 
of NSI, as described in the text. The values of the survival 
probabilities for v e (P SU rv) and i> e (P S urv) for each case are 
also indicated. 



NSI parameters. It is important to pay attention to the 
possible ocurrence of the internal /-resonance and to its 
effect in the external H and L-resonances. The first can 
induce a genuinely new observable effect, item 2 above. 

Here we concentrate on neutral current-type non- 
standard interactions, hence there will be not effect in the 
main reaction in water Chcrenkov and scintillator detec- 
tors, namely the inverse beta decay, D e +p — > e + + n 98 1. 
For definiteness we take NSI with d (down) quarks, in 
which case the NSI effects will be confined to the neu- 
trino evolution inside the SN and the Earth, through the 
vector component of the interaction. 

From all possible combinations of NSI parameters we 
will concentrate on those for which the internal / tran- 
sition does take place, namely |e 7 | > 10 -2 , see Fig. [2] 
Concerning the FC NSI parameters we will consider \e' eT \ 
between few x 10 -5 and 10 -2 , range in which the I- 
resonancc is adiabatic, see Fig.|4] In the following discus- 
sion we will focus on the extreme cases defined in Table Q] 
One of the motivations for considering these cases is the 
fact that the resonances involved become either adiabatic 
or strongly non adiabatic, and hence the survival prob- 
abilities in the absence of Earth effects or shock wave 
passage, become energy independent. This assumption 
simplifies the task of relating the obscrvablcs with the 
neutrino schemes. 



A. Shock wave propagation 

During approximately the first two seconds after the 
core bounce, the neutrino survival probabilities are con- 
stant in time and in energy for all cases mentioned in Ta- 
ble Q] Only the Earth effects could introduce an energy 
dependence. However, at t ~ 2 s the iJ-resonance layer is 
reached by the outgoing shock wave, see Fig. Q] The way 
the shock wave passage affects the neutrino propagation 
strongly depends on the neutrino mixing scenario. In the 
absence of NSI cases A and C will not show any evidence 
of shock wave propagation in the observed D e spectrum, 
cither because there is no resonance in the antincutrino 
channel as in scenario A, or because the //-resonance is 
always strongly non-adiabatic as in scenario C. How- 
ever, in scenario B, the sudden change in density breaks 
the adiabaticity of the resonance, leading to a time and 
energy dependence of the electron antincutrino survival 
probability P surv (E,t). In the upper panel of Fig. [8] we 
show P smv (E,t) in the particular case that two shock 
waves are present, one forward and a reverse one [oo| . 
The presence of the shocks results in the appearance of 
bumps in survival probability at those energies for which 
the resonance region is passed by the shock waves. All 
these structures move in time towards higher energies, as 
the shock waves reach regions with lower density, leading 
to observable consequences in the v e spectrum. 

We now turn to the case where NSI are present, which 
opens the possibility of internal resonances. When such 
/-resonance is adiabatic the situation will be similar to 
the case without NSI. For normal mass hierarchy, AI and 
Cla, v e will not feel the //-resonance and therefore the 
adiabaticity-breaking effect will not basically alter their 
propagation. In contrast, for inverted mass hierarchy 
and large $13, case /?/, the //-resonance occurs in the 
antincutrino channel and therefore v e will feel the shock 
wave passage. However, in contrast to case B now v e will 
reach the //-resonance in a different matter eigenstate: 
D" 1 instead of z?™, sec Fig. [3l That means that before 
the shock wave reaches the //-resonance the v e survival 
probability will be P S urv ~ cos 2 $12 ~ 0.7. Once the 
adiabaticity of the //-resonance is broken by the shock 
wave then v e will partly leave as v™ and therefore the 
survival probability will decrease. As a consequence one 
expects a pattern in time and energy for the survival 
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FIG. 8: Survival probability P SU rv{E,t) for P e as function 
of energy at different times averaged in energies with the en- 
ergy resolution of Super-Kamiokande; for the profile shown in 
Fig. [T] Upper panel: case B is assumed for sin 2 $13 = 10 -2 . 
Bottom panel: case BI, with e TT = 0.07, e eT = 10~ 4 and the 
rest of NSI parameters put to zero. 



probability in the case BI to be roughly opposite than in 
the case B, see bottom panel of Fig. [8l The position of 
the peaks and dips en each panel do not exactly coincide 
as the value of e TT roughly shifts the position of the H- 
resonance. 

In the left panels of Fig. [9] we represent in light-shaded 
(yellow) the range of e eT and s TT for which this opposite 
shock wave imprint would be observable. In the upper 
panels we have assumed a minimum value of the electron 
fraction of 0.06, based on the numerical profiles at t = 
2 s of Fig. Q] In the bottom panels Y~ e min is set to 0.01, 
inspired in the profiles at t = 15.7 s. It can be seen how as 
time goes on the range of e TT 's for which the /-resonance 
takes place widens towards to smaller and smaller values. 
This is a direct consequence of the steady dclcptonization 
of the inner layers. 

For smaller $13, case Clb, the situation is different. 
Except for relatively large e er values the //-resonance will 
be strongly non-adiabatic, as in case C . Therefore the 
passage of the shock waves will not significantly change 
the v e survival probability and will not lead to any ob- 



servable effect. In the right panels of Fig. [9] we show 
the same as in the left panels but for sin 2 $13 = 10~ 7 . 
Whereas for large values of $13, left panels, the H- 
resonancc is always adiabatic and one has only to ensure 
the adiabaticity of the /-resonance, for smaller values of 
$13 the adiabaticity of the //-resonance strongly depends 
on the values of e eT and discussed in Sec. IIVBI 

This can be seen as a significant reduction of the yel- 
low area. Only large values of cither e eT or e Tr would 
still allow for a clear identification of the opposite shock 
wave effects. In dark-shaded (cyan) we show the region 
of parameters for which Ph lies in the transition region 
between adiabatic and strongly non-adiabatic, and there- 
fore could still lead to some effect. 

A useful observable to detect effects of the shock prop- 
agation is the average of the measured positron energies, 
(E e ), produced in inverse beta decays. In Fig. [TUl we 
show (E e ) together with the one sigma errors expected for 
a Megaton water Cherenkov detector and a SN at 10 kpc 
distance, with a time binning of 0.5 s, for different neu- 
trino schemes: case B and case BI with different values of 
e TT . For the neutrino fluxes we assumed the parametriza- 
tion given by Refs. [m,\m with {E (v e )) = 15 MeV and 



(Z?o(^V,t)) = 18 MeV and the followin g ra tio of the total 
neutrino fluxes $o(^e)/ < £'o(^i,T) = 0.8 [99j. 

One can see how the features of the average positron 
energy are a direct consequence of the shape of the sur- 
vival probability, where dips have to be translated into 
bumps and vice-versa. 

Thus, it is important to stress that whereas in case 
B one expects the presence of one or two dips (depend- 
ing on the structure of the shock wave, see Ref [60j), or 
nothing in the other cases, one or two bumps are ex- 
pected in case BI, as seen in the upper left panel of 
Fig. [TO] As discussed in Ref. 60( the details of the 
dips/bump will depend on the exact shape of the neu- 
trino fluxes, but as long as general reasonable assump- 
tions like (Ep c ) < (Ev^ T ) are considered the dips /bumps 
should be observed. 



B. Time variation of Y e 

We have just seen how the distorsion of the density 
profile due to the shock wave passage through the outer 
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FIG. 9: Range of e TT and e eT for which the effect of the shock 
wave will be observed. In the upper panels a minimum value 
of y e mm = 0.06 based on the numerical profiles at t = 2 s 
has been assumed, see Fig. [1] In the lower panels we have 
considered a case with Y e mm = 0.01 inspired in the profile 
at t = 15.7 s. The value of sin 2 #13 has been assumed to 
be 10 -2 and 10 -7 in the left and right panels, respectively. 
We have also superimposed isocontours of constant hopping 
probability 0.1 (blue) and 0.9 (red) in the / (solid lines) and 
H (dashed lines) resonances for inverted mass hierarchy and 
E = 10 MeV and antineutrinos. The area in yellow represents 
the parameter space where both resonances will be adiabatic. 
In the cyan area the /-resonance is assumed to be adiabatic 
whereas H lies in the transition region. 

SN envelope can induce a time-dependent modulation in 
the v e spectrum in cases B and BI. However the time 
dependence of the electron fraction Y e can also reveal the 
presence of NSI leaving a clear imprint in the observed 
v e spectrum, as we now explain. 

As discussed in Sec. IIV Al the region of NSI parame- 
ters leading to /-resonance is basically determined by the 
minimum and maximum values of the electron fraction, 
Y" e min and Yg max . The crucial point is that as the delep- 
tonization of the proto-neutron star goes on, the value of 
y" e min steadily decreases with time. As a result, the range 
of NSI strengths for which the /-resonance takes place 



FIG. 10: The average energy of vp — > ne events binned in 
time for case B (dashed blue) and BI (solid red). In each 
panel different values of e TT have been assumed. The error 
bars represent 1 a errors in any bin. e eT = 10 -4 . 

increases with time, as can be seen in Fig. [2] 

Let us first discuss the observational consequences of 
the time dependence of the electron fraction in case BI. 
If e TT (e 1 in general) is large enough the /-resonance will 
take place right after the core bounce. In this case, as 
seen in the upper left panel of Fig. \TU\ the two bumps we 
have just discussed in Sec. IV Al would be clearly observed. 
However for smaller NSI parameter values it could hap- 
pen that the /-resonance occurs only after several sec- 
onds. In particular for the specific Y e profile considered 
we show how this delay could be of roughly 2, 4 or 9 sec 
for values of e TT of 0.025, 0.02 or 0.015, respectively, see 
last three panels Fig. [101 As can be inferred from the 
figure this delay effect can lead to misidentification of 
the pure NSI effect. So, for instance, in the upper right 
panel, one sees how the two bumps might also be inter- 
preted as two dips, given the astrophysical uncertainties. 
This subtle degeneracy can only be solved by extra in- 
formation on, for example, the time dependence of the 
spectra or the velocity of the shock wave. Given the su- 
pernova model, however, the time structure of the signal 
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could eventually not only point out the presence of NSI 
but even potentially indicate a range of NSI parameters. 

Let us now turn to the normal mass hierarchy sce- 
nario (cases AI and CI a). In analogy to the BI case, 
if e 1 is relatively large the onset of the /-resonance will 
take place early on. As can be inferred from Fig. [3] that 
implies that v e will escape the SN as v 2 . For smaller 
values, though, it may happen that the /-resonance be- 
comes effective only after a few seconds. This means 
that during the first seconds of the neutrino signal v e 
would leave the star as v\ (cases A and C). Then, 
after some point, the electron fraction would be low 
enough to switch on the /-resonance, and consequently 
v e would enter the Earth as T>i. This would result in 
a transition in the electron antineutrino survival prob- 
ability from P surv ~ cos 2 $12 = 0.7 to sin 2 7?i 2 = 0.3. 
Given the expected hierarchy in the average neutrino en- 
ergies (Ep e ) < (E p T ), it follows that the change in Y e 
would lead to a hardening of the observed positron spec- 
trum. The effect is quantified in Fig. [TTJ for different 
values of e TT . The figure shows the average energy of 
the Dp — ► ne + events for the case of a Megaton water 
Cherenkov detector exactly as in Fig. [TQl but for scenar- 
ios AI and Cla. One can see how for e TT = 0.07 the I- 
resonance condition is always fulfilled and therefore there 
is no time dependence. However for smaller values one 
can see a rise at a certain moment which depends on the 
magnitude of e TT . A similar effect would occur in case 
Clb. 



C. Earth matter effects 

Before the shock wave reaches the //-resonance layer 
the dependence of the neutrino survival probability in the 
cases we are considering, on the neutrino energy E is very 
weak. However, if neutrinos cross the Earth before reach- 
ing the detector, the conversion probabilities may become 
energy-dependent, inducing modulations in the neutrino 
energy spectrum. These modulations may be observed 
in the form of local peaks and valleys in the spectrum of 
the event rate c/jP plotted as a function of 1/E. These 
modulations arise in the antineutrino channel only when 
v e leave the SN as i>\ or v 2 . In the absence of NSI this 
happens in cases A and C , where v e leave the star as v\. 
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FIG. 11: The average energy of up — > ne + events binned in 
time for case AI and Cla and different values of e TT . The 
error bars represent 1 a errors in any bin. e eT = 10 -4 . 

In the presence of NSI v e will arrive at the Earth as V\ 
in cases BI, and as v 2 in case AI and Cla. Therefore 
its observation would exclude cases B and Clb. This 
distortion in the spectra could be measured by compar- 
ing the neutrino signal at two or more different detectors 
such that the neutrinos travel different distances through 
the Earth before reaching them 0,0- However these 
Earth matter effects can be also identified in a single de- 
tector 
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By analyzing the power spectrum of the detected neu- 
trino events one can identify the presence of peaks located 
at the frequencies characterizing the modulation. These 
do not dependend on the primary neutrino spectra, and 
can be determined to a good accuracy from the knowl- 
edge of the solar oscillation parameters, the Earth matter 
density, and the position of the SN in the sky 76j. The 
latter can be determined with sufficient precision even if 
the SN is optically obscured using the pointing capability 
of water Cherenkov neutrino detectors 



81] 



This method turns out to be powerful in detecting the 
modulations in the spectra due to Earth matter effects, 
and thus in ruling out cases B and Clb. However, the po- 
sition of the peaks does not depend on how v e enters the 
Earth, as v\ or v 2 . Hence it is not useful to discriminate 
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case AI and Cla from the cases A, C, and BI. 

The time dependence of Y e , however, can transform 
case B into BI, and C with inverse hierarchy into Clb, 
leading respectively to an appearance and disappearance 
of these Earth matter effects. In case BI the presence of 
the shock wave modulation can spoil a clear identification 
of the Earth matter effects. Nevertheless, the disappear- 
ance of the Earth matter effects in the transition from 
case C to Clb allows us to pin down case Clb. 



D. Neutronization burst 

The prompt neutronization burst takes place during 
the first ~ 25 ms after the core bounce with a typical 
full width half maximum of 5-7 ms and a peak luminos- 
ity of 3.3-3.5xl0 53 ergs _1 . The striking similarity of the 
neutrino emission characteristics despite the variability 
in the properties of the pre-collapse cores is caused by 
a regulation mechanism between electron number frac- 
tion and target abundances for electron capture. This 
effectively establishes similar electron fractions in the in- 
ner core during collapse, leading to a convergence of the 
structure of the central part of the collapsing cores, with 
only small differences in the evolution of different pro- 
genitors until shock breakout 
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Taking into account that the SN will be likely to be 
obscured by dust and a good estimation of the distance 
will not be possible, the time structure of the detected 
neutrino signal should be used as signature for the neu- 
tronization burst. In Ref. 78[ it was shown that such a 



time structure can be in principle cleanly seen in the 
case of a Megaton water Chcrenkov detector. It was 
also shown how the time evolution of the signal depends 
strongly on the neutrino mixing scheme. In the absence 
of NSI the v e peak could be observed provided that the 
v e survival probability P VeVe is not zero. As can be seen 
in Tabic U this happens for cases B and C . However 
for case A (normal mass hierarchy and "large" $13), v e 
leaves the SN as v%. This leads to a survival probability 
Pv B v e ~ sin 2 $13 < 10 _1 , and therefore the peak remains 
hidden. 

Let us now consider the situation where NSI are 
prensent. For normal mass hierarchy v e , which is born 
as 14™ passes through three different resonances, /, H 



and L. Whereas / and L will be adiabatic, the fate of 
H will depend on the value of $13. For "large" values, 
case AI, the if-resonance will also be adiabatic. This 
implies that v e 's will leave as > the survival probability 
will be P Ve u e ~ sin 2 $12 ~ 0.3, and therefore the peak 
will be seen, as in cases B and C . If $13 happens to 
be very small, case Cla, then H will be strongly non- 
adiabatic and therefore v e will leave the star as 1/3. As a 
consequence the neutronization peak will not be seen. 

For inverse mass hierarchy, v e is born as v™ and tra- 
verses adiabatically / and L. This implies that they will 
leave the star as v\ and therefore the peak will also be 
observed. However now the survival probability will be 
larger, P VeUe ~ cos 2 $12 ~ 0.7. Thus for a given known 
normalization, i.e. the distance to the SN, one expects a 
larger number of events during the neutronization peak 
in this case. In Fig. [12] we show the expected number 
of events per time bin in a water Cherenkov detector in 
the case of a SN exploding at 10 kpc, for two different 
neutrino schemes, C and BI, and for different SN pro- 
genitor masses. One can see how the difference due to 
the larger survival probability is bigger than the typi- 
cal error bars, associated to the lack of knowledge of the 
progenitor mass. 

Two comments are in order. The neutronization v e 
burst takes place during the first milliseconds, before 
strong dclcptonization takes place. As a result, in con- 
trast to other obscrvables we have considered in this pa- 
per, here the /-resonance will only occur for e 1 > 10 _1 . 
On the other hand in the presence of additional NSI with 
electrons this would significantly affect the v — e cross 
sections, and consequently the results presented here. 



VI. SUMMARY 

We have analyzed the possibility of observing clear sig- 
natures of non-standard neutrino interactions from the 
detection of neutrinos produced in a future galactic su- 
pernova. 

In Sees. IIIII and llVl we have re-considered effect ofv — d 
non-standard interactions on the neutrino propagation 
through the SN envelope within a three-neutrino frame- 
work. In contrast to previous works we have analyzed 
the neutrino evolution in both the more delcptonizcd in- 
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FIG. 12: Number of events from the elastic scattering on elec- 
trons, per time bin in a Megaton water Cherenkov detector for 
a SN at 10 kpc for cases C (dashed lines) and BI (solid lines). 
Different progenitor masses have been assumed: 13 Mq (nl3) 
in red, 15 M (sl5s7b2) in black, and 25 M Q (s25a28) in blue. 
1-sigma errors are also shown for the 15 Mq case. 



ner layers and the outer regions of the SN envelope. We 
have also taken into account the time dependence of the 
SN density and electron fraction profiles. 

First we have found that the small values of the elec- 
tron fraction typical of the former allows for internal NSI- 
induced resonant conversions, in addition to the standard 
MSW-H and MSW-L resonances of the outer envelope. 
These new flavor conversions take place for a relatively 
large range of NSI parameters, namely \e aa \ between 
10~ 2 — 10 _1 , and \e eT \ > few x 10~ 5 , currently allowed 
by experiment. For this range of strengths, in particu- 
lar e rr , non-standard interactions can significantly affect 
the adiabaticity of the H-resonance. On the other hand 
the NSI-induced resonant conversions may also lead to 
the modulation of the v e spectra as a result of the time 
dependence of the electron fraction. 

In Sec. [V] we have studied the possibility of detecting 
NSI effects in a Megaton water Cherenkov detector us- 
ing the modulation effects in the v e spectrum due to (i) 
the passage of shock waves through the SN envelope, (ii) 
the time dependence of the electron fraction and (iii) the 



Earth matter effects; and, finally, through the possible 
dctcctability of the ncutronization u e burst. Note that 
observable (ii) turns out to be complementary to the ob- 
servation of the shock wave passage, (i), and offers the 
possibility to probe NSI effects also for normal hierarchy 
neutrino spectra. 

In Table HT1 we summarize the results obtained for dif- 
ferent neutrino schemes. We have found that observable 
(i) can clearly indicate the existence of NSI in the case 
of inverse mass hierarchy and large #13 (case BI). On 
the other hand, observable (ii) allows for an identification 
of NSI effects in the other cases, normal mass hierarchy 
(cases AI and Cla) and inverse mass hierarchy and small 
$13 (case Clb). Therefore a positive signal of either ob- 
servable (i) or (ii) would establish the existence of NSI. In 
the latter case this would, however, leave a degeneracy 
among cases AI, Cla, and Clb. Such degeneracy can 
be broken with the help of observables (iii) and the ob- 
servation of the neutronization v e burst. The detection 
of Earth matter effects during the whole supernova neu- 
trino signal would rule out case Clb since, as discussed in 
Sec. IV CI a disappearance of Earth matter effects would 
take place due to a transition from C to Clb. Finally, 
the (non) observation of the neutronization burst can be 
used to distinguish between cases AI and Cla. 

Similarly, other degeneracies in Table HT1 may be lifted 
by suitably combining different observables. For exam- 
ple, a negative of observable (ii) could mean either neg- 
ligible NSI strengths or (NU) NSI parameter values so 
large that the internal resonance is always present. In 
this case one could use the observation of the neutron- 
ization burst in order to establish the presence of NSI for 
the case of inverse mass hierarchy. In addition the ob- 
servation of the shock wave imprint in the v e spectrum 
would provide additional information on $13. 

In conclusion, by suitably combining all observables 
one may establish not only the presence of NSI, but also 
the mass hierarchy and probe the magnitude of $13. 
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TABLE II: Expectations for the observables discussed in the 
text: modulation of the P e spectrum due to the shock wave 
passage, the time variation of Y e , the Earth effect, and the 
observation of the u e burst within various neutrino schemes. 
Asterisks indicate that the effect differs from that expected 
in the absence of NSI. See text. 
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